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Abstract 

It has recently been noticed that heavy quark (Q) and gluon (g) fragmentation to 
heavy quarkonium states (QQ) can be calculated in perturbative QCD. This technique 
| allows for the inclusion of a subset of higher order corrections to quarkonium produc- 

tion which dominates at large transverse momentum. A comparison will be presented 
between the new calculation of quarkonium production to the existing data from the 
Tevatron. In addition some new results on predictions for double heavy quark baryons 
at various colliders will be presented. 

in : 

Our understanding of the production and decays of quarkonium states is not perfect. It 
■ was supposed that J/ip production would be dominated by two sources, the charmonium 
production model (CPM)|IJ and B meson decay (BPM)Q. The CDF collaboration mea- 

qq \ sured the J/ip production cross section using their 1988-89 dataH, and found rather poor 
agreement with the sum of the two dominant production mechanisms. They were able to fit 
the data to a sum of these two mechanisms if the normalizations were allowed to float. This 

lO ■ allowed an extraction of the fraction of J/ip mesons from B meson decay, and that, in turn, 
allowed for the inclusion of inclusive J / ip in the measurement of the 6-quark p T distribution. 
The J/tp data points were significantly higher than the theoretical predict ions 0, though 
the UAl collaboration found good agreement between theory and experiment, using similar 
techniques[[5| . This lead to many investigations of the gluon distribution in the proton and 
- the of b- quark pro duct ion ||. 

The questions about the source of the disagreement between theory and experiment were 
settled when the CDF detector was upgraded to include a silicon vertex detector (SVX). 
With the SVX, it was possible to separate the prompt J/ip mesons (i.e., those originating at 
the production vertex) from those produced in B meson decay (i.e., those with a displaced 
vertex). Using their 1992-93 data, the CDF collaboration was able to state that there were 
far too many prompt J/ifc's, and the number produced in B decay was in good agreement 
with theoretical predictions 0. Thus the problem shifted from our understanding of the 
gluon distribution and 6-quark production to that of J/ip production. 

Due to the relatively large mass of the c-quark, ot s is small and heavy quarkonium produc- 
tion and decays can be calculated using perturbative QCD (pQCD), and these calculations 
were, of course, performed. S'-wave state production and decay (J/ip, r\ c , . . .) are well be- 
haved. However, a long-standing problem in quarkonium physics that of the hadronic decay 
(and similarly in the hadronic production) of P-wave states. It was noticed that the matrix 
element of P-wave quarkonium states to 3 gluons has a soft singularity (again, even though 
the S'-wave production and 3 gluon decays are well behaved) ||. Reasonable results are ob- 
tained if one uses the binding energy, confinement radius or the radius of the bound state to 
cut-off the divergences, but these are completely unjustified and non-rigorous procedures. 
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The latter problem was addressed in Ref. || (which was later expanded into a rigorous 
theory called Non-Relativistic QCD, or NRQCDflTO])- A physical hadronic state is made up 
of a superposition of an infinite number of Fock states: 

| M) = ipQQ | QQ) + if) QQg | QQg) + ■■■ (1) 

Generally, the states with one or more gluons don't contribute significantly, as the radiation 
of gluons by heavy quarks Q is suppressed by order v 2 . However, in the case of P-wave decays, 
the angular momentum barrier suppresses the P-wave annihilation by order v 2 , relative to 
the 5-wave annihilation. As the QQ in the | QQg) state can be in the .S-wave, both the 
P-wave color singlet annihilation and the S'-wave color octet annihilation are the same order 
in NRQCD. Thus the decay of the P-wave physical states is given by: 

T(n[ 3 Pj] -> hadrons) = Hi(n)fi(QQf Pj] -> partons) + H 8 (n)T 8 {QQ[ 3 S 1 } -> partons) (2) 

where the subscripts 1 and 8 refer to the color singlet and octet, respectively. The Hi and 
H 8 are non-perturbative parameters. The soft divergence found in earlier works was due to 
the missing H 8 parameter. 
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Fig. 1. Contributions of the various mecha- 
nisms to prompt J/i/j production at the Teva- 
tron. The Xj contributions include the branch- 
ing ratio for radiative decay to J ftp. 



V i 






Direct f Production = 


\ 1 


0; CDF Preliminary 




Dots: LO f I 


■ V\ 1 


Dashes: Fragmentation ifl 


E '■. W 1 


Solid: Total 


E \X s 




- Upper: M=Mo/ 2 




— Lower: /i=2/j, . 




; ^ = V (p T 2 + M 2 ) lor LO 




^ = p T /z for fragm 





5 10 15 20 

Pt (GeV) 

Fig. 2. The total prompt J/ip production rate 
at the Tevatron. The direct (dotted), fragmen- 
tation (dashed) and total (solid) contributions 
are shown separately; also shown are the CDF 
preliminary data points. 



The solution to the J /if) production at CDF problem is due to the existence of pre- 
viously uncalculated fragmentation contributions. These contributions are higher order in 
a s , and so were thought to be negligible compared to the direct production mechanisms 
included in Ref. 0. However, the direct production mechanisms required one to calculate 
QQ production where the quark-anti-quark pair have both small relative momentum and 
the correct quantum numbers to form the relevant bound states. These restrictions modify 
the p T distributions of the produced states. For example, g + g — > g + g falls off as 1/p*, 
while g + g — ► J /if) + g falls off as 1/p 8 , at large p T . The fragmentation of a gluon into 
a quarkonium state will not modify the p T distribution significantly. Including the proper 
factors of a s , the large p T limit of the fragmentation production of J /if) is of order a^/p A 
while direct production goes as a^/pt,] similarly, the large p T limit of \j (the P-wave, spin 
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1 states) production via fragmentation is of order ctf/p 4 while the direct production goes as 
a^/p 6 . At large enough p T , the slower fall-off with p T of the fragmentation contributions can 
make up for the extra powers of a s . The correct calculation of the fragmentation production 
of P-wave states requires the use of NRQCD, and the color octet contribution is important. 
Many fragmentation functions have involving heavy quark states have been calculated, and 



those used in for the following analyses are: gluon fragmentation to quarkonium states pi 
heavy quark fragmentation to quarkonium states[|l!|; photon fragmentation to 
heavy quark fragmentation to heavy quark- heavy quark diquarks[|l4|]. 
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Fig. 3. The total J/ip rate at the Tevatron. 
The prompt (dotted), B meson decay (dashed) 
and total (solid) contributions are shown; also 
shown are the CDF and DO preliminary data 
points. 
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Fig. 5. Pseudorapidity distributions for the 
various contributions to double heavy quark 
baryon production at the Tevatron. 
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Fig. 4. The total ip' rate at the Tevatron. The 
direct (dotted), fragmentation (dashed) and to- 
tal (solid) contributions are shown; also shown 
are the CDF preliminary data points. 
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Fig. 6. Momentum fraction distributions for 
the various contributions to double heavy quark 
baryon production at the Tevatron. 



Now the problem of J /if) production at the Tevatron may be fully addressed. Three 
groups presented compatible results at very nearly the same time|T5|, |ro| |TT| . Results from 
Ref. |TJ] will be presented below. All 2 — > 2 subprocesses were included; MRSDO [18 parton 
distribution functions were used, and // = p T (of the fragmenting parton) for the fragmen- 
tation contributions and /i = M T (of the J /if)) for the direct contributions unless otherwise 
noted (here /x is the fragmentation, factorization and renormalization scale, all chosen to 
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be equal); | r] |< 0.6 was used to simulate the detector acceptance. In Fig. 1, the various 
contributions to prompt J/ip production are shown. In Fig. 2, the total prompt J/ip pro- 
duction rate is shown; the scale /i is varied to show the theoretical uncertainties due to the 
scale choice. The direct contribution is taken from Ref. [19|. The theoretical results are 
compared to the CDF dataf7j. In Fig. 3, the total J/ij) rate is shown, including both prompt 
(direct and fragmentation) and B decay mechanisms; theoretical results are compared to the 
CDF0 and DO 



data. In Fig. 4, the results of ip' are shown, compared to the CDF [21 
data. The total estimated theoretical error (including scale dependence, QCD and relativis- 
tic corrections) is of order a factor of 2. The J/ip production rate is now in good agreement 
with data; the ifr' discrepancy will be discussed later. 
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Fig. 7. p T distributions (at small rapidity) for 
the various contributions to double heavy quark 
baryon production at the Tevatron. 
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Fig. 8. p T distributions (at small rapidity) for 
the various contributions to double heavy quark 
baryon production at the LHC. 



The next set of results to be presented are preliminary, and based upon work in progress. 
Using the fragmentation functions calculated by Falk and collaborators fl4|], the production 
rate for double heavy quark diquarks can be calculated: c — > (cc); b — > (bb); c, b — > cb; and 
c, b — > cb*. The notation is that (QQ) is the spin 1 diquark made up of the same flavor quarks, 
QQ' is the spin diquark made up of different flavor quarks and QQ'* is the spin 1 diquark 
made up of different flavor quarks. It is assumed that the diquark always hadronizes into a 
double heavy quark baryon. The fragmentation probabilities for c — * S cc , S* c ; b — > Ab c ; and 
b — > Efc c , E be* are all of order 10 -5 , while the remaining probabilities are strongly suppressed 
by a factor of (m c /m b ) 3 . The production rates at the NLC (assumed to be a 60 fb~ l /yr, 
500 GeV e + e~ collider) are small, of order 30 events/year. A similar event rate was obtained 
at HERA. Event rates are substantial at both the Tevatron (of order 20k events/year) and 
the LHC (10 7 events/year). In Fig. 5, the pseudorapidity r\ distributions at the Tevatron 
are shown for the various contributions. In Fig. 6, the momentum fraction distributions at 
the Tevatron are shown, demonstrating the strong peaking at large z (z is defined to be the 
ratio of the momentum of the double heavy quark hadron to that of the fragmenting heavy 
quark). Similar results on the r] and z distributions are obtained for the LHC, with the 
primary difference being a smaller relative contribution of the process c — > (cc) at the LHC. 
In Figs. 7 and 8, da / dp T / drj at rj = is shown at the Tevatron and the LHC, respectively. 
Reconstruction of these states will pose a challenge to the experimentalists. Work on this 
project is continuing, and final results will soon be ready. 
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Fig. 9. Different topologies of the lowest- 
order Feynman diagrams contributing to 77 — > 
B r bc. 



Fig. 10. Integrated cross sections for 77 — > 
B c bc versus the cm. energy. The different cal- 
culations are explained in the text. The produc- 
tion mechanisms are as classified in Fig. 9. 



The authors of Ref. |S| calculated the production of B c mesons in photon-photon col- 
lisions, using the Feynman diagrams shown in Fig. 9 (the set (I) can be subdivided into 
(lb) and (I c ) depending on which flavor quark is connected to the photons). They find that 
the recombination diagrams (set (II)) contribute strongly, and that in this case, the heavy 
quark state production can not be well described by a fragmentation process. Fig. 10 shows 
a comparison of the cross section for the full calculation, and for sets (If,) and (I c ) sepa- 
rately, as a function of center of mass energy. Furthermore, the authors of Ref. P3| compare 
the results of set (I) to the relevant heavy quark fragmentation result; the c fragmentation 
contribution does not agree well with set (I c ) while the b fragmentation contribution is in 
relatively good agreement with set (lb). The inclusion of the additional diagrams for the 
gluon-gluon fusion into B c , as well as the convolution of the parton level subprocess with 
parton distributions may change the conclusions at a hadron collider. It should be noted 
that the conclusions of Ref. [^] apply to heavy quark fragmentation contributions only, 
and so may affect the production rate of double heavy quark baryons; the results on J/ijj 
production via fragmentation are dominated by gluon fragmentation, and the argument of 
Ref. |23[] do not apply. After this presentation, another paper discussing the production of 
B c mesons in hadron colliders appeared [24J. The parton level cross sections (g + g — » B c + X) 
are still found to be significantly different, depending on whether the full calculation or just 
the fragmentation approximation is used. After convoluting the parton level cross sections 
with the gluon distributions in the proton, the fragmentation approximation is found to be 
in good agreement with the full calculation for p T > 10 GeV. It is thus likely that the results 
presented here on double heavy quark baryon production at the Tevatron and at the LHC 
will be in quite good agreement with a full (though much more complicated) calculation. 

There remains a large discrepancy between the theoretical predictions for ip' production 
and the CDF experimental data (see Fig. 4). As can be seen from Fig. 1, J/ifj production 
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is dominated by x, T production followed by the radiative decay Xj ~^ J/i> + 7- The ip' is 
more massive than the Xj family, and it is expected that the radial excitation x'j family 
will be above open charm threshold. However, the fall-off of the observed p T distribution 
of ip' appears to be consistent with a fragmentation type production and completely in- 
consistent with a conventional, hard (i.e., direct) production mechanism. Several authors 
suggest a metastable x u i produced dominantly via gluon fragmentation and decaying with 



large branching ratio to ip' as a plausible solution to this problem^, £6], |27j. A search 



strategy for these x'j i n B meson decays at the CLEO II detector has been proposed p8 



Similarly, conventional cc states above open charm threshold, but with quantum numbers 
forbidding the decay into DD have been proposed as a possible solution PB| , as have hybrid 
ccg statesp6[. These states would also be produced dominantly via gluon fragmentation, 
and decay with large branching ratio to ip'. Finally, it has been suggested that color octet cc 
can be produced via fragmentation, and these cc can evolve nonperturbatively to a if)' plus 
hadronsJZg. 

After this presentation, Ref. [[30] appeared. It generalizes the fragmentation mechanism 
and includes additional terms from the NRQCD expansion. Their results are in excellent 
agreement with data for J/ip, ip' and the lowest three T states. 

The author would like to thank NSERC (Natural Sciences and Engineering Research 
Council - Canada) for support, the organizers of MRST at the University of Rochester for 
the opportunity to present this talk and the authors of Ref. p3| for allowing the use of their 
Figs. 1 and 2 (which became Figs. 9 and 10 in this report). 
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